Photoelectron Spectra of Aluminum Cluster Anions: Temperature Effects and 

Ab-initio Simulations 
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Photoelectron (PES) spectra from aluminum cluster anions, Al~ (12 < n < 15), at various tem- 
perature regimes, were studied using ab-initio molecular dynamics simulations and experimentally. 
The calculated PES spectra, obtained via shifting of the simulated electronic densities of states by 
the self-consistently determined values of the asymptotic exchange-correlation potential, agree well 
with the measured ones, allowing reliable structural assignments and theoretical estimation of the 
clusters' temperatures. 
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Photoelectron spectroscopy (PES) is a rich source of 
information pertaining to the electronic structure and 
excitation spectra of atoms, molecules and condensed 
phases. Materials clusters exhibit a high sensitivity of the 
electronic spectrum to the geometrical structure which 
often differs from that of the bulk, and show a high 
propensity to form structural isomers that may dynam- 
ically interconvert at finite temperatures. Consequently, 
high-resolution PES emerges as an important tool in clus- 
ter science, particularly in the face of severe difficulties 
in applying common direct structure-determination tech- 
niques to such systems. 

However, a reliable interpretation of PES spectra is 
often theoretically challenging due to several factors, in- 
cluding: finite-state effects, electronic and ionic relax- 
ations, thermal ionic motions, and structural isomeriza- 
tions. With the advent of accurate ab-initio methods 
for electronic structure calculations, theoretical investi- 
gations of some of these issues have been pursued [Q-^ . 
Particularly pertinent to our study is the development of 
methods which allow practical and reliable simulations 
of PES spectra including dynamical (finite-temperature) 
effects [2-4a]. 

In this paper we address, via the use of ab-initio BO- 
LSD-MD (Born-Oppenheimer local-spin-density molecu- 
lar dynamics) simulations [4a] , methodological issues per- 
taining to simulations and analysis of finite-temperature 
PES spectra. We performed accurate (and practical) cal- 
culations of PES spectra from recorded density of states 
of the clusters using a "generalized Koopmans' theorem" 
(GKT) 1^, concurrent with simulations of the ionic d 
namics. Furthermore, in conjunction with measured 
high-resolution PES spectra for Al~ (12<n<15) clus- 
ter anions, we illustrate that the simulated spectra pro- 
vide a (quantitatively) faithful description of the mea- 
sured ones, including thermal effects, thus allowing re- 
liable assignments of ground as well as isomeric struc- 
tures. Additionally, we demonstrate that through com- 



parisons between simulated PES spectra and those mea- 
sured in three (experimentally undetermined) tempera- 
ture regimes estimates of the clusters' temperatures can 
be obtained. 

In the BO-LSD-MD method the motions of the ions 
evolve in time according to classical equations of motion, 
with the electronic Hellmann-Feynman forces evaluated 
at each MD time step via a self-consistent solution of the 
Kohn-Sham (KS) equations, using the LSD exchange- 
correlations after Ref. Q, and in conjunction with non- 
local norm-conserving pseudopotentials An impor- 
tant element of the method, distingiiishing it from those 
used previously in PES studies (|J^, is the fact that it 
does not employ supercells (periodic replicas of the sys- 
tem), and consequently charged systems as well as those 
having permanent and/or dynamically developed multi- 
pole moments are simulated accurately in a straightfor- 
ward manner on equal footing with neutral ones (i.e. alle- 
viating the need for an artificial neutralizing background, 
large calculational cells, and/or approximate treatment 
of long-range multipole image interactions). 

The ground state structures of Alj~2~-'^lr5i determined 
by us through structural optimization starting from those 
of the corresponding neutral clusters [4b], are displayed 
in Figure 1. Aluminum clusters in this size range favor 
energetically icosahedral-based structures [^; Al^2 
ing an oblate deformed shape, that of Alj~3 being close to 
an ideal icosahedron, and those of Alj~4 and Alj~g being 
capped icosahedra. For Al^g we find that in the energy- 
optimal structure the two capping atoms are located on 
the opposite sides of a " core" icosahedron, resulting in a 
strongly deformed prolate shape (see Fig. 1) pO[ ]. 

The electronic structure of the ground state cluster 
anions exhibits sizable gaps (Eg) between the highest- 
occupied KS molecular orbitals (HOMO) and the lowest 
unoccupied ones (LUMO), as well as odd-even alterna- 
tion (as a function of the number of electrons) in the ver- 
tical detachment energy (vDE) shown in Table 1. Al^g is 
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electronically "magic" (i.e. 40 valence electrons), having 
an exceptionally high vDE |^ and the largest Eg . Its elec- 
tronic structure reflects the corresponding neutral clus- 
ter, which was found [4b] to exhibit a clear jellium-type 
filling sequence ls^lp^l(i^°2s^ for the lowest 20 single- 
particle states. The remaining 20 states, which would 
correspond to jellium l/^^2p^ states, are grouped into 
two broadly overlapping subbands (finite temperature 
broadening of these bands is displayed in Figure 2) and 
show significant pf mixing; this level scheme is known to 
be a consequence of the Ih icosahedral symmetry. 

Although the KS states are not necessarily the "true" 
molecular orbitals of the system, it has been observed 
that the KS HOMO eigenvalue of the TV-electron sys- 
tem, chomo (N) , bears a well-defined relation to the ion- 
ization potential I{N) and electron affinity A{N) JlT|,pt; 
through Janak's theorem jlj]. Following these ideas we 
make use here of a "generalized Koopmans' theorem" 
(GKT) ill] 

eHOMoiN) - C - -IgkHN), (1) 

where is the asymptotic limit of the exchange- 
correlation potential. This nonvanishing energy shift is 
required for an accurate description of the asymptotic 
KS equations ||]. While rigorously the vertical detach- 
ment energy would be given by E{N — 1) — E{N), where 
E{N) and E{N — 1) are, respectively, the total energies of 
the anion and neutral (unrelaxed) clusters, Eq. (1) sug- 
gests a practical approach for evaluation of the threshold 
region of finite-temperature PES spectra through MD 
simulations. Accordingly, neglecting hole-relaxation ef- 
fects, the vDE for removing the electron from the HOMO 
state may be well estimated by {—chomo + v^c)i P^'O" 
vided that remains constant to a good approxima- 
tion, regardless of spatial details of the TV-electron sys- 
tem (such as isomeric atomic configurations of the clus- 
ter). To explore the validity of this condition we have 
calculated for each of the cluster anions the energy shift 
v^^ = E{N-l)-E{N) + eHOMoiN) for a selected set of 
structures (including the ground state one and 10 other 
configurations chosen randomly from finite-temperature 
MD trajectories). The calculated values of were 
found to have a spread of no more than 0.04 eV for each 
of the clusters, and furthermore we found that the de- 
pendence of on N (37 < N < 46) is very weak (see 
Table I). 

While this procedure could be repeated to deter- 
mine values for vDEs from deeper (lower-energy) KS 
states, we have chosen to use a simpler (and more prac- 
tical) procedure whereby we use the shift calculated for 
the HOMO level also for the deeper states In this 
way we generate the full PES spectra from the density 
of states (DOS) recorded in the course of the BO-LSD- 
MD simulations jH]. As shown below, this theoretically 
founded procedure yields a faithful description of the ex- 
perimental data, and it is a viable and reliable alternative 
to previously used methods for finite-temperature ab- 
initio MD simulations of PES spectra which were based 



on either ad-hoc shifts (aligning the theoretical DOS with 
the dominant features in the measured spectra) ]Q or on 
a first-order pertubative treatment [3b]. Furthermore, 
the comparative analysis of the simulated PES spectra 
and the measured ones (see below) validates a posteriori 
certain general assumptions underlying MD simulations 
of PES spectra, i.e.: neglect of finite- lifetime effects of 
the hole (see also Ref. [|| where it is noted that such 
effects may contribute only for very small clusters); the 
use of vDEs (i.e. neglect of ionic relaxations following the 
detachment process); and, assumed equal weights for all 
states contributing to the PES spectrum (i.e. neglect of 
photoelectron transition-matrix effects, which may affect 
line-shape features and certainly the absolute PE cross- 
sections, but not the locations of spectral features i.e. 
binding energies). 

The measured PES spectra for Al]^2~Alr5 ^^"^ shown 
in Figures 2 and 3 (solid line). It is found that clus- 
ters leaving the nozzle early (short residence time) are 
quite " hot" whereas clusters leaving the nozzle late (long 
residence time) are " colder" . Indeed the PES spectra for 
the cold clusters shown in Fig. 2 and the bottom panel of 
Fig. 3 exhibit well-defined features. On the other hand, 
hot clusters exhibit much broader and diffused spectral 
features, as shown in Figure 3 for Al^^g, where we display 
spectra measured for three different residence times, la- 
beled as "cold", "warm", and "hot". Comparisons be- 
tween the locations (binding energies) of the peaks and 
shoulders in the measured and simulated spectra for the 
cold clusters (simulation temperatures of 130 K to 260 K, 
see caption to Fig. 2), validate the u^~shifting procedure 
of the calculated DOS described above. The widths of 
the peaks in the theoretical PES spectra originate solely 
from atomic thermal vibrations since at these low tem- 
peratures isomerization effects and/or strong shape fluc- 
tuations do not occur. The good agreement achieved 
here, without any adjustable parameters other than the 
ionic temperature in the MD simulations, strongly in- 
dicates that the "cold" clusters in the experiments are 
indeed well below room temperature. 

Theoretical PES spectra corresponding to isomeric 
structures of Alj"3-Alj"5, calculated at K, are also shown 
in Fig. 2 (see inset for the threshold regions of Al^g, and 
the dotted line in the panels for Alj~4 and Al^g). The 
isomers for A1J~3 and Alj~4 are the aforementioned deca- 
hedral ones H], and in the Al^g isomer two neighboring 
triangular facets are capped. Comparison between these 
spectra and those calculated for the ground state clusters 
as well as with the measured ones, suggests overall that 
at low temperatures either these isomers do not occur, or 
that their abundance in the cluster beam is rather low. 
In this context we note that starting from the decahedral 
isomer of A1J~3, it transformed readily during short MD 
simulations into the icosahedral one at about room tem- 
perature. This supports our conclusion pertaining to the 
low abundance in the cold beam of clusters "trapped" 
in isomeric structures; however, an even small relative 
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(quenched) concentration of such isomer in the cold Alj~3 
beam may be sufficient to account for the low-binding en- 
ergy tail observed in the measured PES spectra for A\^^ 
(see inset in Fig. 2). 

Both the experimental and theoretical PES spectra, 
shown in Figure 3 for A\^^, which were measured at the 
three temperature regimes mentioned above and sim- 
ulated at the indicated temperatures, exhibit gradual 
broadening and "smearing" of the PES spectral features 
as the temperature increases. We also observe that 
the binding energy of the main peak is rather insensi- 
tive to the thermal conditions, while the line-shape near 
the threshold region (lower binding energies) exhibits a 
rather pronounced temperature dependence. 

The broadening of the spectral features and the (so 
called) " thermal tail effect" near threshold originate from 
the variations of the electronic structure caused by en- 
hanced vibrational motions at the higher temperatures, 
as well as from increased isomerization rates (e.g. in the 
"warm" regime) governed by the free-energy of the clus- 
ter (that is enhanced contributions of lower frequency 
modes to the vibrational entropy p^]), and from dis- 
ordering (" melting" ) of the cluster in the " hot" regime 
(where inspection of the atomic trajectories reveals fre- 
quent transitions between a broad assortment of config- 
urations). Indeed, examination of the vibrational DOS 
of the simulated clusters (obtained via Fourier transfor- 
mation of the atomic velocity autocorrelation functions) 
revealed a marked gradual softening of the clusters at 
the "warm" and "hot" regimes (that is shifting of the 
vibrational spectrum to lower frequencies) coupled with 
increasing overlap between the frequency regions of the 
various modes due to large anharmonicitics. 

In light of the above we judge the overall agreement be- 
tween the simulated and measured spectra and their ther- 
mal evolution as rather satisfactory, and the remaining 
discrepancies (mainly in line-shapes) may be attributed 
to insufficient sampling during the 5 ps MD simulations 
of the thermally-expanded phase-space of the clusters. 

The methodology developed in this study for practical 
calculations of finite-temperature PES spectra, through 
ab-initio MD simulations of aluminum cluster anions 
with no adjustable parameters other than the clusters' 
temperatures, was demonstrated to yield results in agree- 
ment with high-resolution PES spectra measured at var- 
ious thermal conditions of the cluster beam. Such com- 
parative analysis allows reliable structural assignments 
and theoretical estimation of the clusters' temperatures, 
as well as gaining insights into the electronic and struc- 
tural properties of clusters and their thermal evolution. 
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TABLE I. Number of valence electrons (A'^) in the cluster 
anions, HOMO-LUMO gap (Eg), vertical detachment energy 
(vDE) of K ground state anion, and numerically determined 
asymptotic exchange-correlation shift (v^). Energies in eV. 
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FIG. 1. The ground state geometries of Aljj-Aljg (left to 
right). 

FIG. 2. Measured photoelectron spectra of cold (long res- 
idence time) Alj~2- Alj~5 at 193 nm (solid lines) compared to 
the simulated spectra (dashed lines). The simulation temper- 
atures are 160, 260, 200, and 130 K for Alj^^-Aljl;, respectively. 
The arrows correspond to the vDE of the ground state struc- 
ture at K, given in Table I. The inset shows the K PES 
spectra for the ground state (solid line) and the decahedral 
isomer (dashed lino) of Alj^g. The dotted line for Alj4 is the 
K spectrum of the decahedral isomer, and that for Al]~5 is 
the K spectrum of an icosahedral-based isomer (see text). 

FIG. 3. Measured temperature-dependent PES spectra of 
Al]"3 (solid lines) compared to the simulated ones at 930 K, 
570 K and 260 K (dashed lines). HOT - short residence time, 
WARM - medium residence time, COLD - long residence time. 
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